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Abstract 

The supersonic flow held over a body of revolntion 
incident to the free stream is simulated numerically on 
a large, array processor (the CDC Cyber 205). The 
conflguration is composed of a cone-cylinder forebody 
followed by a conical afterbody from which emanates 
a centered, supersonic propulsiTe jet. The free-streara 
Mach number is 2, the Jetrexit Mach number is 2.5, and 
the Jet-to-free-stream static pressure ratio is 3. Both 
the external flow and the exhaust are ideal air at a 
common total temperature. The thin-layer approxima- 
tion to the time-dependent, compressible, Reynolds- 
averaged Navier-Stokes equations are solved using an 
implicit flnite-difference algorithm. The data base, of 
5 million words, is structured in a *pencir format so 
that efficient use of the array processor can be realized. 
The computer code is completely vectorized to take 
advantage of the data structure. Turbulence closure 
is acheived using an empirical algebraic eddyTviscosity 
model. The conflguration and flow conditions cor- 
respond to published experimental tests and the com- 
puted solutions are consistent with the experimental 
data. 

Litrodaetion 

In 1980, a computational study was described in 
which the three-dimensional flow fleld over axisym- 
metric boattailed bodies at moderate angles of attack 
was simulated.^ The exhaust plumes were modeled by 
solid plume simulators, and a second-order-accurate, 
implicit flnite-difference algorithm was used to solve 
the governing partial differential equations on the 
ILLIAC rV array processor. Several flow flelds were 
computed and the results compared with published ex- 
perimental data. The promising resnits of that flrst 
study provided the incentive to extend thr' work to 
include propulsive exhaust jets emanating' rrom the 
afterbody base. The ILLIAC IV was subsequently 
removed from service, however, and it became neces- 
sary to scale down the size and scope of the study to 
the capacity of existing computer resources. 
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In January 1983, the results of a study of super- 
sonic axisymmetric flow over boattails containing a 
centered propulsive jet were presented.^ Those results, 
obtained using a Cray IS computer with 10* words 
of main memory, were compared with existing ex- 
perimental data. Jet-tofree-stream static pressure 
ratio *nH nozzle exit angle were varied parametrically; 
and the predicted trends agreed well with experiment. 

The purpose of this paper is to describe the 
vectorized implementation of the three-dimensional 
Navier-Stokes code on a Cyber 205 computer for boat- 
tailed afterbodies at moderate angles of attack that 
contain a centered propulsive jet. Some computed 
results, which correspond in part to a published ex- 
perimental study for a like conflguration and flow con- 
ditions, are included for illustration. 

Afterbody Conflguration 

The geometric conflguration is a 9 caliber body of 
revolution composed of a 14** half-angle conical nose, 
a cylindrical forebody, and an 3* half-angle conical 
afterbody of 1 caliber length. Centered inside the 
afterbody is a conical nozzle with exit diameter of 0.6 
caliber that is flush with the afterbody base. The 
nozzle exit half- angle is 20**. 

Experimental studies for the same conflguration 
were performed by White and .Agrell® for the model 
immersed in an air stream flowing at = 2.0 and 
a jet-exit Mach number of 2.5. White and Agrell con- 
sidered angles of incidence to the free stream up to 3° 
and jet-to-free-stream static-pressure ratios up to 15. 
Because of limited acces to the Cyber 205 computer, 
computed results are included in this paper only for 
the case in which the angle of incidence is 6** and the 
jet-to-frue-stream pressure ratio is 3.0. 

Governing Equations 

The equations describing the flow are the 
Reynolds-averaged Navier-Stokes equations. These are 
written below in strong conservative form in general- 
ized coordinates as 

dtQ -i- df(F ■ f) + • D + T) = 0 (1) 
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and and?s are the Cartesian unit rectors and 
are the contravariant base rectors, which, 
can be written, aa 

3 * — 4 - tfyZy. 4 - rft'Sa 

3* » SiZ»+ ipiy + &i!«. 

The camponeiita of m om e ntui B^ andpw, 

are in Cartaian. space and the relocitjr rector If 
is generally expressed in terms of the eontravariant 
relocity components, U, V, and W as 

wlr 

-crj^-hn^ + WiTr 

where 3^, 3^, and3f are the corariant base rectors writ> 
ten as 

+ y<?r + 

?i» “ 

3f “t" JTt?^ + r.r2». 

The Jacobian f of the transfonoation is gnes by 

-h ar? y<^r + 


The flax rector F can be decomposed into a 
parabolic part, Fp', which contains only gradient 
diffusire terms, and a hyperbolic part, Fh, which con- 
tains only conrectire-iilce terms, as 


P’h 


f PQ ^ 
puq + pc, 
pvq 4- pcy 
pviq 4- pe» 

V {e + p)q J 


, Fp=^F-Fff (2) 


For flows in which the shear layers are thin (when J?e > 
> 1) and aligned with one principal plane (say the 
plane normal to the q coordinate), the parabolic part 
of F can be neglected in the other two coordinates (f 
and f), without any real loss in accuracy. This is con- 
sistent with boundary-layer theory and yet maintains 
the coupling between the viscous and Lnriscid regions 
that is critical in simulating interactire flows. With 
this thin-layer approximation, Eq. (1) is rewritten as: 


dtQ -h d^iFii; + a^{F • r ) + 9,{Fh • r ) « 0 (3) 

Compntatkinal Grid 

A body-oriented computational grid is constructed 
Insbinaiiiur compatible with the thinrlayer approxima- 
tion. Shown in Fig. 1 is the grid used in the 
present computations. Figure la shows the complete 
eonllgiration and Fig. lb the detail in the base region 
of the afterbody. Radial grid lines on the forebody join 
the surface orthogonaUy. On the afterbody and in the 
exhaust plume^the ra^ai lines are normal to the body 
axis. There are 81 points distributed along the body. 
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Fig. 1 Computational grid: bilateral plane of sym- 
metry. a) Complete configuration (140 x lOO x 20); 
b) Base-region detail. 
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with elostering near the nose and near the base. Of the 
81 points, 21 are used to define the afterbody shape; 
the afterbody is 1 caliber long. An additional 59 points 
are distributed downstream of the afterbody to a dis~ 
tance equal to 21 forebody diameters from the nozzle 
base. These 140 total points define the ^ coordinate 
distribution. The radial distribution, corresponding to 
the Tj coordinate, extends from the body surface to a 
distance equal to 30 forebody diameters both ahead 
of the nose and normal to the body axis. A total of 
60 points is used in this region, with a high degree of 
streching used in order to resolve the sublayer of the 
turbulent boundary layer. (Here the first grid point off 
the body surface corresponds approximately to a value 
of of 8 where An 

additional 40 points are distributed across the nozzle 
and its blunt base, extending from the centerline to 
the body surface. Of these, 20 are in the jet exit plane 
and 20 are on the blunt base itself. 

One- and two-parameter hyperbolic-tangent strech- 
ing- fimctions^ are used in the base region to focus 
resolution near the comers and to achieve a smooth, 
piecewise continuous distribution of points across the 
exhaust plume and base. At the nozzle exit, points are 
distribute along an are describing the conical 
flow exit plane (that is, the arc radius is 
equal to the nozzle exit radius of 0.3 caliber 
divided by the sine of the nozzle-exit half- 
angle of 20°). Downstream of the nozzle, the 
grid lines are aligned so as to clpseiy ap- 
proximate the exhaust plume shape for an ex- 
perimentally observed Asymmetric flow by Agrell 
and White, ^ which is for the same geometric 
configuration and frea-stream conditions, but for a. 
jet-to-free-stream pressure ratio of 9. The 
third dimension, is generated by rotating the 
two-dimensional (^, rj) grid about the cylindrical 
axis while maintaining a uniform angular dis- 
tribution between the rotated planes. Here, 
20 radial planes are used with planes 2 and 
19 coinciding with the bilateral plane of sym- 
metry, where plane 2 corresponds to the lee and 
plane 19 to the windward. Planes 1 and 20 
are image planes used to enforce a symmetry 
boundary condition. Thus, there are (f, rj) 
planes distributed every 10.588° around the half- 
body. 

The total grid dimensions are ( 140 x lOO x 20), cor- 
responding to the f, r/, andf directions, iwspectively, 
for a total of 280,000 points. Of these, (80 x 40 x 20), 
or 64,000, lie inside the body and are not used in the 
computation, leaving an actual total of 216,000 points 
used in the computation. 

Data Structure 

There are 23 variables required at each grid point 
corresponding to the 5 conserved quantities in the Q 


vector, 5 residuals for the solution vector, 9 metric 
coefficients, the Jacobian of the transformation, and 

3 components of vorticity used in the turbulence 
teansport model. This results, for a computational grid 
of 216,000 points, in a data base of 5 x 10* words. 

To accommodate this large data base on a vector 
processor with, a limi ted main memory, the computa- 
tional grid is divided into subsets called ‘blocks.” This 
data structure was originally devised for implemen- 
tation on the ILLIAC IV array processor by Lomax 
and Pulliam and is described in detail in Ref. 6. In 
the present ease, each block is a 20 x 20 x 20 cube 
for a total of 8,000 points and a data base subset of 
184,000 words for the 23 variables. The blocks are 
stacked together in each coordinate direction to form 
a sequence of blocks called ‘pencils.” 

For a given coordinate direction, one complete 
pencU of data is losided into the central memory, and. 
computations are performed on that data correspond- 
ing to the coordinate direction. At any point in the 
computation, only IT variables are required to be in the 
main memory at one time (6 of the 9 metric coefficients 
are not used in any given direction). This results in 
a data-base subset of 136,000 words. For a proces- 
sor with 10* words of main memory then, as many 
as seven blocks of data can be held in storage for im- 
mediate processing. The block dimension is an ad- 
justable parameter and is limited only by the maxi- 
mum pencil length and the main memory of the vector 
processor. 

Shown in Fig. 2, in physical coordinates, are the 
block boundaries for the present configuration. Figure 
2a shows the complete configuration and Fig. 2b the 
detail in the afterbody region. Figure 3 shows the 
corresponding block structure in computational space. 
The mesh nodes of the computational domain are ar- 
ranged in a rectangular latice with positive integer 
coordinates (f, f). Each node belongs to three pen- 

cils, a ^-pencil, an r/-pencil, and a f-pencil. The pencils 
of each sweep direction are given a definite order. For 
the ^-pencils, the 7-coordinate varies most rapidly as 
the pencil index increases; for both the 7-pencils and 
f-pencils the coordinate ^ varies most rapidly. Figure 

4 illustrates this sequencing for the present data struc- 
ture. 

Within a pencil, the planes are naturally ordered 
by the sweep coordinate. The pencils of data can be 
stored in the correct pencil ordering for just one sweep 
direction only. When sweeping in the other direc- 
tions, pencils of data are gathered and fetched for com- 
putation and scattered back when writing the updated 
values. Additionally, the ordering of nodes within a 
plane can be correct for just one sweep direction, and 
it is necessary to transpose the the data in memory 
so that each plane of nodes normal to the sweep direc- 
tion forms a contiguous set of memory locations. In 


189 



the present code, the ordering of nodes Is correct for 
the $>direetioii and transpose rontines are used for the 
other sweep directions. 



{-PCNaL PLANES 



if-PENCIL PLANES 




Fig, 2 Block boundaries; pfaysicai space, a) Complete 
configuration; b) Base-region detaiL 


r-PENOL PLANES 



Fig. 4 Data strnctnre within pencil data base. 
Nninarieal Algorithm 

The numerical algorithm used to solve Eq. (3)' 
ia the approximate factored scheme of Bea'~ and 
Warming.^ Rewriting Eq. (3) as 


n 



Fig. 3 Block boundaries: computational space, com- 
plete configuration. 


d,Q = -di(Ffr-S^)-d^{F -p)-d,{Ffr-^) = Z (4) 

the corresponding difference equation is then 

— (5) 

where the operators are defined by 
L^ = (/ + At<5eA"- 

X., = (/ + - et J-K\i"J) 

Lf =» (/ -f AtSf fl* - €[ y ‘ Vf A^y) 

je« =. ~^6^{jfh ■ tT - Q" 

z^ ^ ~^ts^(jF ■ rT - «E (v,A,)^y g" 
z, == -Ats^iJFH ■ fT - €e y ‘ (Vf Ad="y g" 

where the 5^, 5,^, and 5^ are central-difference 

operators; V,, and Vj are backward-difference 
operators; and A^, A,,, and A^ are forward-difference 
operators in the f-, 7-, and ^-directions, respectively. 
The At term is a forward-difference operator in time. 
For example, 
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Begin Loop 


— <3(f + Af, 7. f) - «(?. 7, s) 

(.pencils: 


and 

Read: 

Q, J, R, u, c-metrics 

v <<3 - Qis. 7, f) - Qi€ - Ai, 7. f ) 

lYanspoie: 

Q. J, R, « 

The Jacobian matrices 

Compute: 

R=-j^« + Je„ 

. 



A » • g ) 

Transpose: 

R, u 

B *= d^iFu • g^) 

Write: 

R, 0 / 

C=^d,,iF„-T) 

7 -pencils: 



M=r^dQ{Fp-T) 

an described in detail by Pailiam and Steger.' Fonitb* 
order explicit terms (preceded by the coefficient eg) 
and second-order implicit terms (preceded by the 
coefficient f /) have bera added to control nonlinear In- 
stabilities. 

Equation (5) is solved in three' successive sweeps 
of the data base, each sweep inverting one of the 
operators on the left-hand side: 


Read: Q, J, R, ui, tf-metrics 

TVanspose: Q, J, R, u 
Compote: « =•= w(€) + w(f) -|- w(j;) 

Mw) 

R— + 

lYanspose: L^^(R) 

Write: t;'(R) 


L^L^^tQ =« L‘^ ( ^e-h>e, ) 

Li AtQ = r/i-V ( ^< + je, + ) 

AtQ = i-Vr/r; ( -h ^ J 

The solution is advanced in time by adding. AtQ to Q 
after the ^ sweep. 

In the general case, pencils of data are loaded into 
central memory fonr tim^and operated on for each 
time-step advance: once each for the f and 7 direc- 
tions and twice for the f direction. First the right- 
hand side of Eq. (5) is formed and then the left- 
hand-side operators are inverted one by one. A flow 
schematic showing the ordering of operations, includ- 
ing data reads, transposes, computations, and data 
writes is shown below where the symbols R and uj 
represent variables used to accumulate the right-hand- 
side elements and vorticity elements, respectively, for 
each coordinate direction. 

^-pencils: (initial step only) 

Read: Q, J, f-metrics 

Compute: R = Z^, w =* 

Write: R, w 


(-.pencils: 

Read: Q, J, I.li‘(R)> f-metrics 

■ftanspose: Q, J , X.‘‘(R) 

Compute: Lj'L^‘(R) 

Transpose: l;‘l;‘(R) 

Write: l;‘x.;‘(R) 

^pencils: 

Read: Q, J, f-metrics 

Compute: AtQ, Q, R = Z^, w = uiO 

Write: Q, R, u 

End Loop 

bi this flow sequence, 62 variables are read, 57 
variables are transposed, and 31 variables are written. 
Fbr the special case in the present study in which the (- 
pencils are just one block long, a more efficient operar 
tion sequence csm be used that substantially reduces 
the number of reads and writes required. This is shown 
below. 

f-pencils: (initial step only) 


Read: Q, J, |-metrics 

Compute: R = Z^, (J = u(^) 
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Be^ Loop 
^.pencils; 


Read: ^-metrics 

Transpose: Q, J, R, u 
Compute: R 

(J s u(0 4- w(f) 
Transpose: R, u 
Wnte: R, u 


^•pencils: 


Read: 

Transpose: 

Compute: 


Q, J, R, u, 7 -metrics 
Q, /, Rr w 

«=«c*<a4-w(ri + w(7) 
ftiiu) 




^.pencils: 


Read: 

TYansposer 

Compate: 

IVanspose: 

Write: 


f-metrics 

Q. J, ilV(R) 


^pencils: 

Read: Q, J, ^-metrics 

Compute: ^tQi Q 

w = w(0 

Write: Q, 

End Loop 

In this flow sequence, 32 variables are read, 52 are 
transposed, and 18 variables are written, a savings of 
nearly 50% in the I/O. In both the general case and the 
special case, the data read-transpose sequence and the 
transpose-write sequence can be replaced by/tjie more 
efflcient “gather" and “scatter" commands available 
for the Cyber 205 (Ref. 9). Further improvements in 
efficiency can be obtained by using asynchronous I/O 
in conjunction with a rotating memory backing store. 
The most efficient code, however, will be realized by 
using a solid-state backing store in conjunction with 
gather and scatter commands or with a code that is 
fully core contained. 

The numerical algorithm conforms well to large 
vectorization. For block sizes of 20 x 20 x 20, the vector 


length is 400. Timing studies with the present code in- 
dict an MFLOP rate (million of floating- point opera- 
tions per second) oni5 when computing in half preci- 
sion (32-bit word lengths) on a 2-pipe configuration. 
On a 4-pipe configuration the MFLOP rate increased 
to 20T. There are approximately 3,800 floating point 
operations executed for every grid node per time step 
resulting in a CPU time of 33 x 10"* sec per point per 
time-step on a 2-pipe machine and 18 x 10"* sec per 
point per time-step on a 4-pipe machine. The transpose 
times (transposes do not contain any floating-point 
operations) are 5.6 x 10"* sec per point. Equivalent 
transposes performed by gather and scatter instruc- 
tions require just 1.8 x 10"* see per point. When 
synchronized I/O to and from rotating backing store 
was used, the average I/O time was 25 msec per vari- 
able per block. This translates directly into 172 x 10^* 
sec pCT point, but overlapping the I/O reduces this to 
94 X lU^ sec per point. (The Cyber 205 used for these 
timing- studies was configured with four I/O channels 
ta accommodate overlapping.) This time, a result in 
large part of the latency time in accessing disk files, 
can be reduced to. nearly zero by using I/O buffers in 
conjunction with asynchronous I/O or with solid-state 
backing storage. The use of I/O buffers, however, im- 
pUea the availability of additional main memory and 
imposes an additional constraint on the pencil size. To 
avoid this constraint, the data flow should be modified 
such that a subset of contiguous blocks of data in a 
pencil are operated on while blocks at each end of the 
subset are being buffered in and out. 

Boundary Conditions 

Boundary conditions are imposed at the ends of 
each data pencil; the data pencils are identified by 
number in Fig. 3. For the ^-direction, pencil No. 1 
starts at the jet-exit plane. Supersonic conical flow 
conditions corresponding to a jet-exit Mach number of 
2.5 and a static pressure of 3poo are imposed at the 
first data plane. At the last plane of each of the five 
f-pencils, which correspond to the outflow boundary, 
first-order extrapolation is used so that d^Q = 0. 
Pencil No. 2 in the ^-direction begins at the blunt base. 
Here slip conditions and an impermeable adiabatic wail 
are imposed so that 

d^p) — d^pv) — d^{pw) = 0 
pu = 0 

d^[e — 0.5(pu^ pv^ + = 0 

Pencils 3, 4, and 5 in the ^-direction begin on the 
grid centerline of revolution (at ^ = 0) ahead of the 
forebody nose. Here a second order extrapolation to 
the centerline is used such that 

<5«(p) = 5?(p«) = d^pw) = d^{pe) = 0 
while the lateral momentum is set to zero 
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pw «0 

In addition, at eacii 17, the Q yaloes are aTeraged over 
f on the centerline and naed aa boundary Talnea for all 
f at each 7. Special treatment of the bate comer at 
the afterbody-blnnb-base jnnetion is osed to acconnt 
for the singniar nature of that line. For the {’-aweepa, 
the f'line of data in peneU No. 3 that eorrespondi to 
this comer is treated in the same manner aa the lirst 
plane of data in pencil No. 2 that corresponds to th» 
blunt base. This Hue of data is treated differently in 
the 7>aweep and is described in the second paragraph 
following. 

After the forebody flow field is fhUy developed 
during the course of the solution, the first two 7-pencils 
can be dropped from the computation and boundary 
conditions imposed on the {^pencils that correspond 
to the fully dereloped flow at the plane that is the 
upstream tenndary o£ 7-peneil No. 3 . This reduces 
the total data base by six blocks without altecing:the 
ralldity of the solution. This simplification is strictly 
Talid only for supersonic external flows. The solution 
downstream can be further developed to steady state, 
and jet parameters can. even be varied to generate ad- 
ditional solutions. 

Boundary conditions for the 7-direction consist of 
the imposition of free stream conditions at the last 
plane of each of the seven 7-pencils; no-slip, adiabatic 
wall condition for the first plane of 7-peneils. 1 through 
4 , which correspond to the body surface; and first-order 
extrapolation to the centerline for pencils 5 , 6, and T 
snch that d^Q —Q~ Centerline averaging, as deseribed 
for the f-pencil bonndary^ahead of the body, is also 
used for the 7-pencil boundary in the jet. The line of 
data in 77- pencil No. 5 , which corresponds to the comer 
between the afterbody and the blunt base, is treated 
in the same manner as the first plane of 7-pencils I 
through 4 . As a result, this line of data is double 
valued: one value for the ^ sweep described previously 
and the noslip, adiabatic value for the 7-sweep. 

For the f-direetion, bilateral symmetry is imposed 
by setting the data at the first and last ^-planes equal 
to the values in the third plane and in the second from 
last plane, respectively, with a sign change included in 
the lateral momentum component (/>v). 

Torbulenc* Closure 


the eoeSleient of molecular viscosity, { jt, » ft + 
Itr), thereby relating turbulent transport directly to 
gndienU of the mean-flow variables. La a Cartesian, 
coordinate system, the three-dimensional molecular 
stress tensor can be written as 

Ti ““ (y + Os)V» + Tsf + TatiaZ, 

TyaZyZa (p 

TfsZsZg - 4 * “f* (P “h <^a)2a?a 

In the thin-shear-layer approximation, the only com- 
ponents of the stress tensor that are retained are those 
having gradients with respect to 7 only. 

Thrbnient heat transport is defined iu terms 
of meaiMnergy gradients and an eddy-condnctivity 
eoefileient K, such that K, ^ K + Kt- Typically, 
the eddy^ondnetivity coefficient is related to the eddy- 
viscosity coefficient via a turbulent PrandtL number 
Pry where 

Ptt “ CfinlKr 

The turbulent Prandtl number is assumed constant at 
a value of 0 . 9 . 

The algebraic eddy-viscosity model used here is 
that proposed by Baldwin and Lomax.^** This model 
is partieniarly well suited to complex flows that con- 
tain regions in which the length scales are not clearly 
defined. It is described briefly as follows: For wall- 
boonded shear layers, a two-layer formulation is used 
such that 

Mr ” G^T)mi»»r for 7 ^ 7er»»»»»«r 

Mr “ fOV 7 ^ 7«r»M»v«r 

where 7 is the normal distance from the wall and 
is the smallest value of 7 at which values 
from the inner and outer formulas are equal. The 
Prandtl- Van Driest formulation is used in the inner (or 
wail) region. 

(MT)iini€r ~ P^\^\ 

I -» 0.47 (1 — «p(— 7/A)1 

A = 2 QHail^/p^ 

The formulation for the outer region is given by 

(MT)*«t«r * 0"0168C'cyPw«*«Pffj«*(7) 


The Reynolds stresses and turbulent heat-finx 
terms have been included in the stress tensor and 
heat- flux vector by using the eddy- viscosity and eddy- 
conductivity concept, whereby the coefficients of vis- 
cosity aind thermal conductivity are the sum of the 
molecular (laminar) part and an eddy (turbulent) part. 
Eddy-viscosity models incorporate turbulent transport 
into the molecular-transport stress tensor by adding 
the scalar eddy-viscosity transport coefficient nr to 


^mas ^ most \ 

7mos ^ifl ^ maaj 

The quantities 7 mo> Fmas are determined from 
the function 

P( 7 ) » 7 (wj [1 — exp(— 7 /A)] 

where F„aa is the maximum value of F{r]), and r]max is 
the value of 7 at which it occurs. The function Pki«6(7) 
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is the Klebanoff intennittescy fnnction gnen by 

'Hie quantity is the dillbrence between the mazi- 
nrant and minimnin total sqoared Telocity in the profile 
(along an f^-coordinate line), 

“ 9mtf» 

and for boundary layers, the minimnm is defined as 
zero. The other constants are grren by 

1 0,25 , ^Kl^k 0,3 

The adrantage of this model for boondary-iayer 
flows are as follows: 1) tor the inner region, the Telocity 
and length, scales are always well defined, and the 
model is consistent with the “law of the wall”; 2) in the 
outer region for weU-behsred (simpie> boundary layers^ 
where there is a weltdefined len^ scale the 

Telocity scale is determined by Fmaat which is a length 
scale times a vorticity scale; 3) in the outer region of 
complex boundary layers where the length a wall 
becomes meaningless, a new length scale is determined 
from a Telocity dfrided by » Telocity gradient 
( [(>;( }, and the veloei^ scale is q^i/. 

The outer formulation, which is independent of q, 
is also used in the free-shear flow regions of separated 
flow and in regions of strong Tiscous/inriscid inter* 
action. In these- regions the ran Driest damping: 
term, [exp( — q/A)}, is neglected. For jets and wakes, 
the IQebanoff intermittency factor is determined by 
measuring f^om the grid centerline, *nd tiie nuninmm 
term in q 4 i/ Is eralnated fh>m^the profile instead of 
being defln^ as zero. 

The vaiidity of the eddy-viscosity model constants 
for high-pressure, compressible exhaust jets has not 
been established, and compressibility effects are not 
accounted for. 

At the exhaust-jet exit plane and in the uear-base 
region, the eddy viscosity is assumed to be negligibly 
smalt and to increase spatially to the value given by 
the outer model over a short distance downstream of 
the base. 

Computed Reiulta 

As mentioned in a preceding section (Afterbody 
Configuration), a flow field has been computed for 
the body placed at an angle of incidence of 6” to 
a free stream at Mach 2. The jet-exit Mach num- 
ber is 2.5 with a static pressure 3 times that of the 
free stream. Beginning with an impulsive start in a 
uniformly flowing stream at Mach 2, the solution was 
advanced timewise to a dimenstoniess time (tdfUoo) of 
5.1, where d is the forebody diameter and is the 
undisturbed free-stream speed. Although a solution 


at a time of 5.1 is probably not sufficiently eonverged 
to pennit valid qaantitative comparisons with experi- 
ment, it is sufficient to establish the basic flow-fleld 
character and to illustrate the features of the solntion 
and the computer code. 

The initial time-step size of Af ^.0001 was in- 
creased to At «0.001 as the solution passed through its 
initial rapid transient. A variable time-step was used 
in the subsonic flow regime downstream of the base in 
order to minimize the growth of nonlinear instabilities 
aggravated by changes in sign of the eigen-values in 
this region. The time-steps in this subsonic region were 
sealed down by a factor equal to the local streamwise 
Mach number with a entoff minimnni factor of O.OOt 
imposed to prevent the time-step from, going to zero. 

Occurring' physically in this region is a rapid ex- 
pansion of the jet around the nozzle lip followed im- 
mediately by a strong- recompression in the form, of a 
birrei shock; in addition there is a slip surface defining 
the bonndnry between the exhaust plume and the ez- 
temai flow. Each of these three high-gradient features 
is focused at the nozzle lip and demands a high degree 
of resolntion that has not been provided for in the com- 
putational grid used here. 

Shown in Fig 5 are computed density contours 
in the bilateral plane of symmetry in the vicinity of 
the body. The lower surface is the wind side. Clearly 
defined downstream of the afterbody is the slip snrr- 
fhce demarcating the boundary between the exhaust 
plume and the external flow. The propulsive jet ex- 
pands rapidly around the nozzle lip and can induce ‘low 
separation on the afterbody surface. For low-pressure 
jets, or no jet at all, there will be a region of recir- 
culating flow on the blunt base. The afterbody drag is 
strongly influenced by the detail of the separated flow. 



Fig. 5 Coraputeddensitycontonrs, plane of symmetry; 
Afco = 2, Af/ = 2.5, PjiP^ = 3, 
a — 6®, Rci = 1.5x10*. 
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Flf: 6 AAorbodx flair dstail: amtaem streamliius aad dnuitT^ eontonn oa bilatanl plan* of STnunatr^'. 


Tlie detail of tbe separatios patters is shows in 
FI|f. 6 in which eompsted surface streamlines have 
been mapped on the afterbodj and projected on the 
bilateral piane-of-symmetry view of the density eonr 
tour plot over the aft portion of the body onij. There 
is a separation node on the iee generator of the eoni- 
eal afterbody at x ^ 3.92. All surface streamlines on 
the iee side of the body flow into this node. A line of 
separation extends from this node, downward on the 
afterbody surface, to a separation saddle at*s « 3.98, . 
33* rtom the wind generator. The flow direction along 
this line of separation is upward flrom the saddle to the 
node. There is also flow optward (him the separation 
saddle downward to the end of the base, around to the 
wind generator. 

Shown in Fig. T is a perspective view of the 
surface streamlines on the afterbody and the blunt 
base. The outer edge of the base is a diTidlng surface 
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Fig. 7 Perspective view of surface streamlines over 
conical afterbody and annular base. 


streamline extending fhim a saddle point on the lee 
generator to a node point approximately 33* from the 
wind generator. A dividing streamline .can be seen cir* 
enmscribing the annular base connecting a saddle point 
on the- windward and a nodal point on the Iee. This 
line separates the external flow frL'in the flow from the 
jet. Flow is npward ffom the windward saddle to the 
lee-side node. 

Shown in Fig. 3 is a sketch of an end-view projec- 
tion of the full view of the afterbody (not to scale) 
showing all the dividing streamlines and their cor- 
responding singuiar points and flow directions. 



Fig. 3 End-view schematic of dividing surface and 
singular points streamlines. 
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Tbe trajectories of the fluid particles in the piaae 
of symmetiT in the base region are shown in Fig. 9. On 
the lee, seen in Fig. 9a, the fluid lri>m the jet expands 
around the noxzle lip and mores outward toward, the 
edge of the base. Upon, meeting the external flow, it 
tarns downstream deflnes the exhaast plame boai^ 
dary. A region of rererse flow can be clearly seen abore 
the afterbody lee generator. The path of the flaid in 
the external flow is orer this separation region and 
aroand the afterbody base to the slip surface deflning^ 
the boundary between the exhaust plume and external 
flow. The point defined by the outer edge of the base 
and the afterbody lee generator is a singular point that 
from the fluid streamlines, appears as a saddle point in 
both the circumferential plane and in the radial plane, 
and as a nodal point in the streamwise bilateral plane 
of symmetry (tte plane of the base). 

On the windward, shown in Fig. 9b, the stream^ 
lines just off the wind generator of the afterbody 
turn the comer and move toward the slip sntfhee 
between^ the jet. and the external Sow. Ail external 
flow stremnlines (exeiadbig the surface streamline) ap- 
proach the slip surface downstream of a saddle point 
in the bilater^ plane of symmetry located at x « 
9.016 on the plnme-extemal flow boundary. The sni^ 
face streamline turns the comer and approaches the 
windward saddle point on the base itseff. Fluid from 
the jet expands around the nozxle lip and moves out- 
ward The fluid just off the lip mores to the saddle 
point on the base and the flnid farther inside the lip 
expands toward the plume boundary downstream of 
the saddle point on the slip surface. 



Surface-pressure distributions orer the afterbody 
surface and over the base are shown in Figs. 10a 
and. 10b, respeetirely. An expansion at the forebody- 
afterbo(fy junction over the afterbody surface can be 
seen, lliis expansion is greatest on the windward, 
where the pressure level is highest, and decreases 
toward the lee. The cirenmferential variation of pres- 
sure near the lee side is quite small for the entire length 
of the afterbody. Toward the end of the afterbody 
there is a slight recompression on the lee side which 
is not obeerved on the windward. Just at the end of 
the afterbody there is an expansion as the flow turns 
around the afterbody toward the base. 

Figure 10b shows a projected view of the base 
and jet-exit pressure distribution. The left side of the 
*top hat” pressure distribution corresponds to the lee, 
aad the far side corresponds to the windward. The 
large uniform pressure distribution of the top hat” 
configuration corresponds to the high-pressure jet, and 
the undnlating' *brim” of the hat is the distribution on 
the mitwiIm base. On the windward there is a rapid 
expansibn at the nonle Up followed by a fairly large 
recompression toward the outer edge of the base. The 
same trend is observed at other radial positions around 
the hue but to a lesser de g ree. The circumferential 
variation of base pressure is consistent with the ex- 
perimentally ob s e r ve d variation of White and Agreil 
for the same jet-to-ftee-stream pressure ratio. It is 
interesting to note, however, that in most experimen- 
tal studies the radial variation of pressure is assumed 
negligible and is not measured. The distribution in Fig. 
lOb ekariy indicates a substantial variation across the 
aanolar baae^ 



Fig. 9 Base-region path lines: plane of symmetry, a) Lee; b) Windward 
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FIs: la Sufacft pmsnre dlstribotioai perspeetirB. s) Conieai afterbodT; b) Aaanjar baie and jet exit plans. 


Cooeindinc Remark* 

As implicit solntios p ro ced u re’ for the thin> 
lajer appraximatioii to the three-dlmenaional, time> 
dependent, compressible, Reynolda- a Te rag e d Nsrier> 
Stokes equations on a large amj processor has been 
described. An example problem ivas simnlated on the 
Cjber 20S eompnter that required a data- base of 5 x 
10* -words. The eflleient treatment of this- large data 
base has been described is some detail. 

The flow-fleld simnlated was the supersonic flow 
orer a bodT- of rerolntios at incidence to the ftee 
stream. A propuisire jet emanated from the boattailed 
afterbody, inducing a complex, three>dimenaional 
separated-flow pattern. This separated flow-fleld, 
which contributes substantially to the afterbody drag, 
has been described in detail for the particnlar geometry 
and flow conditions considered. The computed solu- 
tion is consistent with experimental data observed for 
the same configuration and flow conditions. 
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